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Dependent  Server  Vacations 
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ABSTRACT 

This  paper  shows  how  to  determine  the  stationary  distribution  of  the 
virtual  wait  in  M/G/1  queues  with  either  one  at-a-time  or  exhaustive  server 
vacations.  Each  vacation  time  may  depend  on  the  immediately  preceding 
service  time  or  on  whether  the  server  finds  the  system  empty  after  return¬ 
ing  from  vacation.  In  this  way,  it  is  possible  to  model  situations  such  as 
long  service  times  followed  by  short  vacations,  and  vice  versa.  The  method 
of  analysis  employs  level-crossing  theory.  Detailed  examples  are  given  for 
various  cases  of  service  and  vacation-time  distributions  ^  . 
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INTRODUCTION 


Recent  interest  in  queues  with  server  vacations  has  focused  on  stocheistic  de¬ 
composition  properties  for  the  number-in-the-system  random  variable.  Re¬ 
searchers  who  have  investigated  and  obtained  a  number  of  important  results 
for  this  aspect  of  the  model  include  Gaver  (1962),  Cooper  (1970),  Levy  and 
Yechiali  (1975),  Scholl  and  Kleinrock  (1983),  Fuhrmann  (1984),  Fuhrmann 
and  Cooper  (1985  a,b),  Shanthikumar  (1988),  and  Harris  and  Marchal  (1988). 
Models  have  also  been  developed  for  G/G/1  server  vacation  queues  in  Doshi 
(1985),  Keilson  and  Servi  (1986),  and  Servi  (1986).  But  few  results  are  avail¬ 
able  for  the  stationary  distribution  of  the  waiting-time  random  variable  in 
such  queues.  The  mean  and  variance  of  the  waiting  time  in  the  specialized 
case  of  D/G/1  queues  with  server  vacation  were  given  in  Servi  (1986)  using 
complex  transform  methods.  Doshi  (1986)  showed  an  integral  equation  for 
the  probability  density  function  (pdf)  of  the  stationary  virtual  wait  for  an 
M/G/1  server-vacation  model,  based  on  Brill  and  Posner  (1977,  1981).  But 
his  equation  did  not  involve  the  kind  of  dependency  considered  in  this  paper, 
and  Doshi  also  did  not  provide  any  of  the  detailed  complexities  necessary  for 
carrying  out  such  analyses. 

A  model  with  state-dependent  server  vacations  was  investigated  by  Har¬ 
ris  and  Marchal  (1988).  Their  paper  offered  a  decomposition  principle  for 
M/G/1  queues  in  which  the  server- vacation  times  depend  on  the  number  of 
customers  remaining  in  the  system  at  service  completion  epochs. 

The  object  of  the  present  paper  is  to  derive  the  stationary  distribution  of 
the  virtual  waiting  time  in  M/G/1  queues  in  which  each  server  vacation  time 
may  depend  on  the  immediately  preceding  service  time.  Consideration  of 
this  dependency  allows,  for  example,  the  possibility  of  rewarding  short  service 
times  with  long  vacations,  penalizing  long  service  times  with  short  vacations, 
or  modeling  situations  where  the  server  naturally  requires  a  service-time- 
dependent  recovery  period  after  each  service.  The  modeling  technique  and 
solution  method  used  in  this  paper  are  based  on  level  crossing  theory  for 
obtaining  probability  distributions  directly,  originated  by  Brill  (1975),  and 
further  elucidated,  for  example,  in  Brill  and  Posner  (1977,  1981). 

The  level-crossing  methodology  is  extremely  advantageous  in  the  present 
context,  and  yields  a  general  integral  equation  for  the  pdf  of  the  virtual  wait 
immediately  (given  in  Section  2,  Equations  (1)  -  (4)).  This  paper  emphasizes 
how  the  general  integral  equation  can  be  specialized  in  various  important 
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cases  of  service-time  and  server-vacation-time  distributions.  It  will  become 
clear  (Section  3)  that  this  procedure  may  require  careful  analysis,  accord¬ 
ing  to  the  type  of  dependency  of  server  vacation  on  customer  service  time. 
Furthermore,  we  obtain  explicitly,  or  indicate  how  to  obtain,  the  pdf  of  the 
virtual  wait  in  these  cases.  The  examples  presented  point  to  a  potentially 
large  class  of  models  that  can  be  solved  in  a  similar  manner. 

Section  2  of  this  paper  derives  a  general  model  equation  for  the  stationary 
distribution  of  the  virtual  wait  in  M/G/1  queues  with  one-at-a-time  server 
vacations  depending  on  the  immediately  preceding  service  times.  Variations 
of  this  equation  are  given  for  the  exhaustive  server  vacation  model  and  other 
special  cases.  Section  3  utilizes  these  model  equations  to  obtain  the  sta¬ 
tionary  distribution  of  the  virtual  wait  in  four  examples.  These  illustrative 
examples  demonstrate  an  approach  for  obtaining  the  distribution  of  the  vir¬ 
tual  wait  in  other  special  cases  of  interest. 

2  Model  Equations  for  the  Stationary  Dis¬ 
tribution  of  Virtual  Wait  in  M/G/1  Queues 
with  One- at- a- Time  Server  Vacations  De¬ 
pending  on  Service  Times 

Customers  arrive  in  a  Poisson  stream  at  rate  A  and  are  served  in  order  of  ar¬ 
rival  by  a  single  server  who  may  take  intermittent  vacations  following  service 
completions.  The  service  times  are  denoted  hy  Sn,  n  >  1,  having  common 
probability  distribution  function  B{  ).  Immediately  following  the  nth  service 
completion,  the  server  takes  a  vacation  of  duration  =  V„,  which  may 

depend  on  S„-  Denote  an  arbitrary  steady-state  service  time  by  5  *=*  5n, 
and  a  corresponding  vacation  time  by  V  V„,n  >1.  If  the  server  re¬ 
turns  from  a  vacation  and  finds  the  system  empty  (no  waiting  customers),  it 
immediately  takes  another  vacation  whose  duration  is  denoted  by  Vq. 

We  define  the  conditional  distribution  function  of  the  random  variable  V 
given  S  by 

Hix,y)^PiV  <x\S  =  y), 
and  the  distribution  function  of  Vo  by 

Ho{x)  =  P{Vo  <  x). 
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Also  define  B{x)  =  1  —  B{x),H{x)  =  1  —  H{x),  and  Ho{x)  =  1  —  Ho{x),  for 
all  X. 

Let  Vr(t)  denote  the  virtual  wait  at  time  t,  defined  as  the  delay  of  a 
potential  arrival  at  time  t  to  the  system  until  he  begins  service  (at  instant 
+  The  stationary  distribution  function  and  density  function  of  W{t) 
are  denoted  by 


F{x)  =  ^lim  PiW{t)  <  x) 

and 

= 

respectively,  assuming  they  exist.  In  this  system  all  customers  must  wait 
a  positive  time  with  probability  one,  and  thus  P(W(t)  =  0)  =  0,  <  >  0, 
implying  that  liTnt^ooPi^V{t)  =  0)  =  0.  Although  general  conditions  for 
ergodicity  are  not  examined  in  this  paper,  these  conditions  are  given  for  the 
special  cases  illustrated  below  in  Section  3.  We  now  employ  system  point- 
level  crossing  theory  to  derive  the  model  equations  for  the  pdf  f{w),w  >  0. 

Level- crossing  theory  uses  sample-path  properties,  and  a  sample  path  of 
the  {Lr(t)}  process  is  shown  in  Figure  1.  Without  loss  of  generality,  assume 
the  system  starts  empty  and  the  server  begins  a  vacation  at  time  t  =  0.  The 
sample  path’s  slope  is  -1  for  all  t  >  0,  except  at  arrival  epochs  and  the  starts 
of  empty  system  server  vacations,  which  are  epochs  of  jump  discontinuity. 
At  each  arrival  epoch,  a  jump  is  generated  equal  to  the  sum  of  the  arrival’s 
service  time  and  the  immediately  subsequent  server- vacation  time. 

From  system  point-level  crossing  theory,  the  long-run  expected  down¬ 
crossing  rate  of  level  ty  >  0  in  the  state  space  is  equal  to  f{w)  (Brill  1975,  or 
Brill  &  Posner  1977,  1981).  The  long-run  expected  upcrossing  rate  of  level 
w  is  equal  to  the  expression  (for  w  >  0) 

X  r  P{S^V  >w-  z)f{z)dz  f{0^)P{Vo  >  w) 

J  *=o 

=  A (/;/„•  B'(-»  -  j  -  y.y)dB{y)\  f{z)dz 
+  >.JZ„B[w-z)J{z)dz  +  f{0*)H„{w), 
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recalling  that  f(0'^)  equals  the  long-run  expected  rate  of  level  zero  hits  (Brill, 
1975). 

Equating  the  long-run  expected  rates  of  sample-path  down-  and  up- 
crossings  of  state-space  level  w  >  0  yields  the  model  equation  for  the  sta¬ 
tionary  pdf  f{w),  w  >  0: 


tVJ  — «  _ 

f{w)  =  a/  /  H(w  -  z -y,y)f(z)dB(y)dz 

J  t=0  •/y=0 

fw  _  _ 

+  X  B(w- z)f(z)dz  + (1) 
J  z=0 

The  normalizing  condition  is  given  by 

/  f{w)dw  =  1  (2) 

Jo 

since  the  virtual  wait  is  concentrated  on  (0,  oo)  with  probability  one. 

If  the  vacation  times  {Kv}  are  independent  of  the  service  times  {^n}  for 
all  n  >  1,  then  Equation  (1)  becomes  (for  w  >  0) 


f(w)  =  xj  I  H{w  -  z  ~  y)f{z)dB{y)dz 

J *=0  Jy—0 

X  r  B{w-z)f{z)dz-\-f{0^)Ho{w),  (3) 

J  *=0 

where  H{x)  =  P{V  >  x),  x  >  0. 

In  the  exhaustive  server-vacation  problem,  vacations  occur  only  when 
the  server  finds  no  customers  in  the  system,  so  that  =  0,  n  >  1,  and 
Equation  (1)  simplifies  for  w  >  0  to 

fW  _  _ 

f(w)  =  X  B{w  -  z)f{z)dz  -f  (4) 

J  z—0 
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Examples 

This  section  contains  four  examples  which  illustrate  how  to  obtain  the  sta¬ 
tionary  distribution  of  the  virtual  wait.  In  the  first  example,  service  times 
have  a  generalized  hyperexponential  {GH)  distribution  (Botta,  Harris 
Marchal,  1987),  and  server-vacation  times  are  exponentially  distributed  with 
means  depending  on  the  duration  of  the  immediately  preceding  service  times. 
In  the  second  example,  the  service  times  are  exponentially  distributed  and 
the  vacation  times  are  deterministic  with  value  depending  on  the  duration 
of  the  immediately  preceding  service  times.  Example  3  is  a  special  case  of 
Example  1  in  which  service  times  and  server-vacation  times  are  assumed 
independent.  Finally,  Example  4  illustrates  an  exhaustive  server- vacation 
model  which  is  a  special  case  of  Example  2. 

Example  1 


Assume  that  each  service  time  is  distributed  as  random  variable  5  having 
the  GH  distribution  with  probability  density 

=  Ec.e-'“-  (x>0), 

t=i 


dP{S  <  x) 
dx 


where  ^,  >  0  and  — oo  <  c^  <  -foo  for  t  =  l,...,n  (n  >  1).  Next  assume 
that  the  server  vacation  times  are  exponential  random  variables  with  means 
1/7(5),  where  7(5)  is  a  function  of  the  immediately  preceding  service  time 
S,  and  there  exists  some  real  number  A  >  0  such  that 


7(5)  = 


( 


71 

72 


if  S<  A 
if  5  >  A. 


In  this  model,  A  may  be  considered  to  be  a  system  control  parameter. 

Let  Vo  be  an  exponential  r.v.  with  mean  1/70-  Then  model  Equation  (1) 
for  the  stationary  pdf  of  the  virtual  wait  becomes 
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f\ii  rw-z 

/(^)  =  aE/  / 

“  Jt=0  Jy=0 


,—yi(w-t-y)  -ti,y 


c,  e  f{z)  dy  dz 


(0  <  w  <  A), 


(5) 


t=i 


and 

^  f\u-A  r  rA 

/(«^)  =  >E/  / 

^  y2=o  [4=0 


rti;  — i 

dy  + 

Jy=A 


^--rj(tu-*-y)  -/i,v 


c,e  '^’*'dy 


f{z)dz 


f{z)dz 


(6) 


-foty 


^  /tu  r  fw^x 

— /  L»/y=:0 

+-'E  /"  (<=./('.)e-"'<”->/(2)*  +  /(0*) 

t=l  •^*=° 

(/I  <  w). 


The  pdf  f{w)  is  continuous  for  all  tt;  >  0  (Brill,  1987)  and  thus  f{A  )  = 
Equations  (5)  and  (6)  can  be  rewritten  respectively  as 


.=1 


LM.(7i  -  M.) 


f{w)  =  A71  ^  ' 

n 


Jt=0 


t=i 


c. 


7i  -  M. 


f 

^«=C 


(0  <  w  <  >1) 


(7) 


and 
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f{w) 


1=1 

+^E 

t=l 

n 

-AE 


7i  - 


Jt=0 


72 


_  r“'“ 

—  /I,  y*=o 


t=l 


72  -  M 


>*.)1  rw-A 

,  yz=o 


4-Ae[-^1  r  -e-^‘ /(2)d2 

[71  -/^.J  ^  ^ 

^  /‘U» 

+  A  V  -  /  e-''‘(“'-*)/(2)  rf2  +  /(0+ ) 

^  M,  ^*=0 


(A  <  w).  (8) 


We  can  solve  for  /(to)  in  (2),  (7)  and  (8)  by  converting  (7)  to  a  differential 
equation,  substituting  its  solution  into  (8),  and  applying  (2).  Alternatively, 
we  may  use  numerical  methods  or  level-crossing  estimation  (Brill,  1987). 

Example  2 


In  this  next  example,  there  is  a  control  parameter  A  >  0,  and  the  vacation 
times  are  deterministic,  depending  on  the  preceding  service  times  or  whether 
the  vacation  originates  when  the  system  is  empty.  The  service  times,  {S„,  n  > 
1},  are  iid  exponential  random  variables  with  mean  l//x.  Let  S  denote  an 
arbitrary  service  time.  If  5  <  A,  the  server  takes  a  vacation  of  duration 
V  =  £>2,  and  if  5  >  A,  the  server  gets  a  vacation  of  duration  Di  <  D2.  Thus 
long  service  times  (5  >  A)  are  followed  by  short  vacations,  and  vice  versa. 
Moreover,  assume  that  vacations  when  the  system  is  empty  are  Vq  =  Do, 
and  that  0  <  Do  <  Di  <  <  A  holds.  Other  cases  for  the  value  of  A  can 

be  treated  similarly.  The  model  equations  for  /(■),  the  stationary  pdf  of  the 
virtual  wait,  are  as  follows,  noting  that  the  function  f{w)  has  a  discontinuity 
at  ID  =  Do  (Brill,  1987): 
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(0  <  ti;  <  Do) 


(9) 


f{w)  =  x[  f{z)dz  +  f(0'^) 

Jz:=0 

f{w)  =  A  /  f{z)dz  {Do  <w<  D2),  (10) 

Jx=0 

f{w)  =  \r  f{z)dz-\- X  r~^\->^^'^~^^~^'>f{z)dz  (11) 

J  —  J  z~0 

(Z>2  <  w  <  A  +  Di), 

f-MJ  rw  —  D^  _ 

f{w)  =  X  /  f{z)  dz-^X  dz 

J  z=.w—D^  Jz^w—A  —  Di 

A  —  Di  -  . 

^  fi^2)dz  (12) 


{A  +  Di  <  tu  <  +  D2), 


fW  fW-Dj 

f{w)  =  X  f{z)dz  +  A  /  e-'^(“'-^'-*>/(z)<f2 

J  i=\u~Di  J  t=vi—A-D\ 


+  A  r  (13) 

Jt-xu  —  A-D’}  ^  ^ 


rw-A-Di 

+  A  /  e-^^'^-^'-‘^f{z)dz  {AaD2<w). 

J  x=0 


We  can  solve  for  f{w)  in  (2)  and  (9)  -  (13)  by  converting  to  differential 
equations,  substituting  the  solutions  for  lower  values  of  w  into  the  equations 
for  higher  values,  and  finally  applying  (2).  Thus,  from  (9)  and  (10),  we  get 


f{w) 


ae^'"  (0  <  lu  <  Do) 
,  6e""'  {Do<D2) 
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where  a  and  6  are  constants  to  be  determined.  Balancing  down-  and  up- 
crossing  rates  at  state-space  level  Do  yields  f{Do)  =  A  J^o  f{z)  dz,  result¬ 
ing  in  6  =  a(l  —  This  procedure  can  be  continued  in  Equations  (10) 

-  (13),  using  the  functional  form  of  f{w)  for  0  <  tu  <  D2,  etc.  Alterna¬ 
tively,  we  may  use  numerical  methods  or  level-ciossing  estimation  to  obtain 
f{w),w  >  0. 

Example  3 


Consider  a  situation  where  the  vacation  times  {Vn,n  >  0}  are  iid  ex¬ 
ponential  random  variables  with  mean  1/r/,  and  the  service  times  are  iid 
exponential  random  variables  with  mean  1/fi,  fi  ^  r}.  The  vacation  and 
service  times  are  independent,  so  that  Equation  (3)  applies;  thus  the  model 
equation  for  /(•)  is 

fW  fW—t 

Jxi=:0  Jy^O 

+  X  r  4-/(0+ )e-’’'^  {w  >  0), 

Jx^O 


which  simplifies  to 


f{w) 


r  e~^{'—)dz  -  r  f{z)dz 

—  ii  Jz=o  “n  —  u  Jt=o 


-H/(0+)e-’''" 


{w  >  0). 


(14) 


Equation  (14)  is  a  special  case  of  (7)  and  (8)  in  Example  1,  with  n  =  1,/ii  = 
A'.Ti  =  72  =  letting  A  =  0  in  (8)  or  letting  A  — ♦  00  in  (7). 

Taking  derivatives  with  respect  to  w  in  (14)  and  rearranging  terms  yields 
the  linear,  second-order  differential  equation 

f"{w)  -f  (77  -f  p  -  X)f'{w)  -f  {rjfi  -  At;  -  A/i)/(w)  =  0.  (15) 
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We  assume  that  the  stationary  distribution  of  the  virtual  wait  exists,  and 
hence  the  mean  interarrival  time  exceeds  the  mean  jump  size  for  non-zero 
waits,  as  in  the  classical  M/G/1  queue,  that  is,  when 

1  1  1 

-r  >  h 
X  7} 

It  follows  that  the  characteristic  equation  of  (15)  has  two  real  roots  with 
positive  product  and  negative  sum,  and  hence  both  roots  are  negative.  De¬ 
noting  the  roots  by  and  R2  yields  the  solution  for  the  pdf  of  the  virtual 
delay  as 


f{w)  =  +  be^^'"  {w  >  0), 


where  a  and  b  are  constants.  Substituting  this  solution  into  (2)  and  equating 
the  coefficient  of  e~^'^  on  the  right-hand  side  of  identity  (14)  to  zero  yields 


with  solution 


_a - 1  _j>_-  =  0 

M+Ri  ^  t*+R7 


I*{Ri-R7) 
m(R7-Ri) 

Example  4 


This  example  treats  a  case  of  the  exhaustive  server- vacation  M/G/1 
model,  so  that  =  0,  n  >  1.  Assume  that  server- vacation  times  when  the 
system  is  empty  are  deterministic,  each  equal  to  Do,  i  e.,  V“o  =  Do-  Moreover, 
assume  for  the  sake  of  exposition  that  service  times  are  also  deterministic, 
with  each  service  S  =  Di  and  Dt  >  Do-  (See  Figure  2.) 
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Equation  (1)  now  becomes 


rw 

f{w)  =  a/  f{z)dz  {w  >  Di) 

f{‘w)  =  A  /  f{z)dz  (Do  <  W  <  Dj) 

J  f=.0 

f^W 

f{w)  =  A/  f{z)dz  +  f(0'^)  (0  <  ty  <  Do). 

Jt=0 

For  w  >  Di,  assume  that  f(w)  =  ce^'“ .  Then 

ce^"'  =  a|  [e^”'  -  , 

so  that 

Ae-^^>  +  ^  -  A  =  0,  (16) 

which  yields  the  value  of  0,  for  <  0,  by  the  uniqueness  of  the  solution  of 
the  integral  equation. 

For  w  €  (Do,  Di), 

f{v})  =  A/(«;), 

so  that 

/(to)  =  de^”'  (d  =  constant). 

By  continuity  at  Di  (Brill,  1987),  /(Dj)  =  ce^^‘,  and  it  follows  that 

d  =  e(^-^>^‘c.  (17) 

For  to  <  Do, 

/'(to)  =  A/(to), 

so  that  /(to)  =  ae^”'  and  /(O'*')  =  a  (a  =  constant).  Note  that  there  is 
a  discontinuity  &t  w  =  Dq  (Brill,  1987).  However,  /(Dq  )  is  the  expected 
entrance  rate  to  [0,  Do);  the  expected  exit  rate  from  [0,  Do)  is 


/•Do 

A  /  f{z)dz 

Jt=Q 


and  therefore 


f{Di)  .=  X  f(z)dz. 
Jx=0 


Thus 


and 


/Do 

=  A  /  ae^^dz, 

Jt=0 


a  ~  c 


exp[{l3  —  \)Di  +  \Do] 


exp{XDo)  —  1 


=  cA. 


(18) 


Finally,  from  the  normalizing  condition  (2),  we  get 

c  r  e^'^dw  +  c  [e(^-^>^*]  e^'^dw  +  cA  e^'^dw  =  1 

Jw^Di  ^  ^  Jw—Do  Jw^O 


or 


g0Oi  e^Do  _  j 

ir-^  A 


=  1.  (19) 


The  solution  is  obtained  by  substituting  the  values  of  c,  a,  d  and  /0  given 
in  (19), (18),  (17)  and  (16),  respectively,  into  the  formulas  for  the  pdf  /(w) 
over  the  state  space  intervals  (0,  Do),(iAo,  iAi)  and  {Di,oo). 
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Figure  3.  PDF  of  virtual  wait  in  the  exhaustive  server  vacation  model. 


4  CONCLUDING  REMARKS 

Level-crossing  theory  is  particularly  useful  for  modeling  and  analyzing  server- 
vacation  queues,  since  service  times  and  vacations  are  additive  in  the  virtual 
waiting-time  process.  Moreover,  the  theory  is  concerned  with  balancing  rates 
of  down-  and  up-crossings,  regardless  of  their  qualitative  origins. 

Examples  1  and  2  (Section  3)  indicate  possible  new  ways  to  control  such 
queues  by  treating  the  parameter  .4  as  a  controllable.  Control  could  be  ex¬ 
tended  by  using  a  collection  of  (positive  integer)  M  such  parameters  {A^}, 
where  0  <  .i4i  <  412  <  ■■■Am-  Clearly,  the  examples  of  Section  3  demonstrate 
solution  techniques  that  can  be  applied  to  a  large  class  of  M/G/l  server- 
vacation  models,  having  various  distributions  of  service  and  vacation  times, 
with  disciplines  possibly  other  than  first-come,  first-served. 
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